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Paths are important structural elements in complex networks because they are finite (unlike 
walks), related to effective node coverage (minimum spanning trees), and can be understood as 
being dual to star connectivity. This article introduces the concept of random path applies it for 
the investigation of structural properties of complex networks and as the means to estimate the 
longest path. Random paths are obtained by selecting one of the network nodes at random and 
performing a random self-avoiding walk (here called path-walk) until its termination. It is shown 
that the distribution of random paths are markedly different for diverse complex network models 
(i.e. Erdos-Renyi, Barabasi- Albert, Watts-Strogatz, a geographical model, as well as two recently 
introduced path-based network types), with the BA structures yielding the shortest random walks, 
while the longest paths are produced by WS networks. Random paths are also explored as the means 
to estimate the longest paths (i.e. several random paths are obtained and the longest taken). The 
convergence to the longest path and its properties ire characterized with respect to several networks 
models. Several results are reported and discussed, including the markedly distinct lengths of the 
longest paths obtained for the different network models. 

PACS numbers: 89.75.Fb, 02.10.Ox, 89. 75. Da 



'The longest journey starts with a single step. ' ( Old 
Chinese proverb.) 

I. INTRODUCTION 

A good deal of the attention so far focused on com- 
plex networks (e.g. i, i, i, i) has been motivated 
by properties involving high degree nodes (hubs in scale 
free models) and shortest paths between pairs of nodes 
(small world networks). The former model is character- 
ized by scale free degree distribution, which accounts for 
the heterogeneity of the local connectivity of nodes and 
the presence of hubs. Scale free models are also small- 
world because each hub establishes shortcuts among all 
nodes to which it is connected to. Recently [1, 0], the 
concept of a star was suggested in order to represent the 
network motifs defined by the edges connected to each 
node. As discussed in those works, a star can be under- 
stood as the dual of a path 16]. Particularly long paths 
can therefore be related to hubs through such a duality. 
Stars and paths, inter-related by the path-star transfor- 
mation and its inverse [6, 7], present completely opposite 
intrinsic properties: stars are centralized, while paths are 
sequential. 

While great attention has been given to the degree (e.g. 
presence of hubs) and shortest paths, other types of paths 
— such as random and longest — have received little at- 
tention. It is particularly interesting to study paths be- 
cause: (i) they have finite length (unlike walks); (ii) they 
correspond to the connected subgraph with the smallest 
possible number of edges linking a set of nodes (a concept 
related to that of minimum spanning tree); and (iii) as 
outlined above, they can be understood as duals of stars, 
a fundamentally important motif in complex networks. 

Random paths, as suggested in this work, can be ob- 



tained by selecting a node at random and performing 
a random path-walk (i.e. a walk not allowed to repeat 
nodes and edges) until unvisited edges or nodes are no 
longer available. Such paths are therefore closely re- 
lated to the exploration Q, sampling [l3| and cov- 
erage [ll| of the nodes in a network while trying to 
minimize the number of edges crossed (not a single edge 
or node in a path are visited more than once). Actually, 
it is interesting to observe that maximal paths, such as 
Hamiltonian paths [l3| covering all the nodes of a net- 
work, are special cases of minimum spanning trees [l^. 
Random paths can also be understood as being produced 
by moving agents which destroy the edges and nodes as 
they are visited, as it could happen during the spreading 
of a disease or attack on a network. In addition to being 
related to such dynamical processes, paths can be con- 
sidered in order to get insights about the structure and 
topological properties of networks. For instance, if the 
several random paths obtained from a given network are 
found to have short lengths, it may be inferred that the 
network contains many bottlenecks (e.g. hubs or edges 
with large betweeness centrality) . Therefore, the analysis 
of random paths can provide additional valuable informa- 
tion about the structure of the analyzed networks. 

The longest path in a complex network corresponds to 
the path (or one of the equal-length paths) involving the 
largest number of edges (and therefore nodes). Alter- 
natively, the longest path can be defined as the random 
path involving the largest number of edges. While the 
distribution of the shortest path lengths bears important 
implications for communications and distributions tak- 
ing place in the network, the longest path is naturally 
related to effective coverage of the nodes by a moving 
agent. For instance, the presence of a Hamiltonian path 
indicates that all nodes of the network can be visited 
while crossing the minimal number of edges. Therefore, 
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the existence of long paths in a specific network may indi- 
cate that these structures are intrinsically important for 
the dynamics or that the network has been constructed 
around or by incorporating such paths (e.g. [tJ). One of 
the possible reasons why scant attention has been focused 
on longest paths is that, while shortest paths can be de- 
termined in polynomial time, longest path identification 
is an NP-complete problem (e.g. |13j). 

The current article aims at investigating random and 
longest paths in diverse theoretical complex networks 
models. It starts by describing the basic concepts and 
methodology, including a simple stochastic method for 
estimating the longest path, and proceeds with the ex- 
perimental investigation of the random and longest paths 
in six theoretical complex networks models. 



II. BASIC CONCEPTS 



A undirected complex undirected network, formed 
by N nodes and E edges among those nodes, can be 
fully represented in terms of its adjacency matrix K 
of dimension N x N. Each existing edge implies 
K{i,j) = K{3,i) = 1, with K{i,j) = K{j,i) = 
being otherwise imposed. Two edges are said to be 
adjacent whenever they share one of their extremities. 
A walk corresponds to any sequence of adjacent edges 
(ii, 22); (^2, is); • ■ • (*p-i, where p is the length of the 
walk. A walk which does not repeat any edge or node is a 
path. An Hamiltonian path [h^ is a path that encompass 
all nodes in the networks. The length of a path is equal 
to the number of its constituent edges. The shortest path 
between two nodes is defined as one of the paths between 
those nodes which has the smallest length. 

The immediate neighbors of a node i are those nodes 
which are connected to i through shortest paths of length 
1. A star is a motif containing a node i and the edges at- 
tached to it. The degree of a node is equal to the number 
of edges in its respective star. The node degree averaged 
within a network is called average degree. The clustering 
coefficient of a node i is the ratio between the number 
of undirected edges between the immediate neighbors of 
i and the maximum possible number of undirected edges 
among those nodes. 

Six theoretical models of complex networks are consid- 
ered in the present work including four traditional mod- 
els — Erdos-Renyi (ER), Barabasi- Albert (BA), Watts- 
Strogatz (WS) and a geographical model (GG) — as well 
as two recently introduced knitted types of complex net- 
works Q — the path-transformed BA model (PA) and 
path-regular networks (FN). The ER, BA and WS net- 
works are grown as traditionally done (e.g. [l!, [3, i, i, H ) . 
The GG network is obtained by randomly distributing N 
nodes within a square and connecting those nodes which 
are closer than a parameter d. The PA and FN net- 
works are grown as described in Q. Basically, the PA 
networks (path-transformed BA networks) are obtained 



by star-path transforming all nodes in an original BA 
network. The FN (path regular networks) model, on the 
other hand, is easily obtained by defining paths involving 
all network nodes in random order and without repeti- 
tion. 

All networks have N = 100 and to = 3 or m = 5 (m 
is the number of spokes in the added nodes in the BA 
model), which implies average degree (k) = 2m. Because 
the average degrees considered in this work are above the 
percolation critical density, most of the nodes in each net- 
work belong to the largest connected component, which 
is the only part of the network considered for the analyses 
reported in this article. 



III. METHODOLOGY 



This section describes the simple stochastic algorithm 
for estimation of the longest path used this article as 
well as the several path-based measurements which can 
be applied in order to characterize the global connectivity 
of networks. Given a network, the longest path estima- 
tion approach involves finding many random paths while 
keeping track of the longest obtained path length. De- 
spite the fact that the determination of the longest path 
is NF-complete, the method is verified, at least for the 
considered networks, to converge steady and quickly to 
the longest path. 

Given a complex network, one may choose a node i at 
random and start a self-avoiding random walk (or path- 
walk [1^) from two outbound edges of i, until unvisited 
edges and nodes are no longer available to the walker at 
both sides of the walk (a situation which will be called 
here saturation or termination of the path). It is pro- 
posed in this work that such a random path can provide 
valuable information for complex networks studies. More 
specifically, provided several random paths are obtained 
by using such a simple algorithm, the distribution of their 
properties, especially their length, can provide valuable 
information about the global structure of the network 
connectivity. 

The algorithm adopted for picking a random path is 
simple and involves the following steps: 

(i) Select a node i at random; 

(ii) Chose two edges A and B emanating from i (in case 
that node has degree one, chose only one edge); 

(iii) Perform a path- walk along both edges A and B until 
the walk stops (i.e. a node is reached which leads to no 
unvisited node or edge). 

The length of a random path is henceforth represented 
by p. 

The algorithm proposed here for the estimation of the 
longest path is also very simple and involves obtaining 
several random paths and taking the one exhibiting the 
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largest length as the longest path. Given that the deter- 
mination of the longest path is NP-complete, it is impor- 
tant to consider the convergence of such a simple algo- 
rithm. This is performed experimentally in Section [IV Bl 

IV. RESULTS AND DISCUSSION 

The experimental investigations reported in this sec- 
tion, reported in the following subsections, involved 50 
realizations of each model (PA, FN, ER, BA, WE and 
GG) with N = 100 nodes and considering and m = 3 
and m = 5. 



A. Random Paths in Network Models 

The first important point to be investigated regards 
the properties of random paths extracted from the sev- 
eral theoretical models of complex networks considered 
here, (i.e. PA, PN, ER, BA, WS and GG). In order to 
understand random paths in these several models, ran- 
dom paths were extracted from a network (iV = 100 and 
m = 3) of each type by using the algorithm described 
in Section IIIII and their lengths p determined. Figure [T] 
shows the length of the first 200 obtained random paths. 
Although this figure considers just one randomly chosen 
network of each type, it was verified that similar results 
are obtained for other network samples, except for the 
WS network, which tended to produce diverse signatures 
for each specific network. It is clear from this figure that 
the lengths of the random paths varied substantially for 
each type of network. The first interesting result is that 
the values of p present relatively small variation, except 
mainly for the FA and ER models; the BA, FN and GG 
structures tended to produce random paths with similar 
lengths. The GG and BA networks led to the shortest 
random paths, while the ER and PA were characterized 
by medium size lengths. The longest lengths were ob- 
tained for the PN, FA and WS models. Observe also 
that the standard deviation of the random path lengths 
tend to increase with the respective average values. 

Before interpreting the variations of the random path 
lengths for each model, an additional, more systematic 
study was performed considering the mean and stan- 
dard deviation of the length p of 2500 random paths for 
N = 100 and m = 3 and m = 5. Figure [2] shows the re- 
spectively obtained scatterplots. Observe that each point 
in these scatterplots was obtained from the measurement 
of the mean and standard deviation of p considering 2500 
random paths from each of the 50 realizations performed 
for each configuration (i.e. N = 100 and to = 3, 5). 

A series of interesting effects can be inferred from these 
scatterplots. Let us discuss first the case m = 3, shown in 
Figure [UJ a). The BA networks, shown in filled diamonds, 
yielded the smallest average value of p, followed by the 
GG, ER, FA, PN and WS models. The BA structures are 
also characterized by relatively small standard deviations 



of p. In other words, paths randomly extracted from BA 
networks tend to be short (about 30 edges long) and with 
similar lengths. On the other hand, random paths in WS 
networks tend to be long (from 70 to 100 edges) and with 
widely varying lengths. The PA and ER networks pre- 
sented remarkably similar measurements, indicating that 
these two models are very similar as far as random path 
properties are concerned for m = 3. The PN structures 
resulted remarkably uniform, in the sense that the mean 
and standard deviation of p resulted particularly similar 
for each of the 50 realizations, defining a compact cluster 
in the scatterplot in Figure [D The GG networks yielded 
widely varying random path lengths extending from 10 
to 80 edges. 

We now proceed to the scatterplot obtained for to = 5, 
i.e. (k) = 10. Such an increase of average degree im- 
plied some changes relatively to the situation discussed 
for TO = 3. First, the random path lengths increased for 
all network models. The BA cluster, for instance, shifted 
from values in the range [15,45] to [25,65]. The largest 
increase of random path lengths was achieved for the GG 
model, which now extends to the maximum length of 100 
and over. Observe that, because of the way in which 
they are grown, GG networks may sometimes exceed 100 
nodes (see Section |lT]). Another interesting effect caused 
by the average degree increase regards the separation of 
the FA and ER cases, with the respective clusters now 
presenting substantially less overlap. Also, the disper- 
sions of the FN and WS networks were substantially re- 
duced. 

Let us know try to discuss each of the main findings 
above, which is done as follows: 

(a.) Longest random paths are obtained for larger av- 
erage degree: Usually, the higher the average number 
of edges in a network (reflected in a high average node 
degree), the more possibilities the agent performing the 
path-walk will have to keep proceeding to unvisited nodes 
and edges. Observe that the extreme case where the net- 
work is fully connected (i.e. each node is linked to all 
other nodes in the network), all random paths will al- 
ways have maximum length of N . 

(h) Short random paths in the BA model: Consider 
a typical BA network, with its respective hubs. Once 
started, a path-walk is very likely to pass through hubs. 
Once such hubs are visited, they can no longer be part 
of a path. Because several of the network edges are con- 
nected to hubs, the paths soon become saturated. 

(c)Long random paths in the WS model: Recall that 
the WS model is grown by starting with a regular network 
(a ring) where each node is connected to a number of pre- 
vious and subsequent nodes. Therefore, by distributing 
the connectivity almost regularly amongst the nodes, the 
chances of having at least one unvisited edge leading to 
an unvisited node increases substantially. Interestingly, 
being small world does not imply short random paths. 
On the other hand, the large variances of p obtained for 
this network model are accounted by the rewirings used 
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FIG. 1: Lengths of the first 200 randomly selected paths obtained for each of the considered networks models. 



to induce the small- world effeet. Because several of the 
paths in this model encompass all nodes, this type of 
network can be understood to frequently contain Hamil- 
tonian paths, implying high efficiency of node coverage 
by moving agents. 

(d) Uniformity of random path lengths in the PN 
model: As shown in [2], PN networks are highly regu- 
lar regarding several topological measurements. Such an 
enhanced topological regularity tends to produce simi- 
lar random path lengths. Still, the substantially distinct 
values of average and standard deviation of the random 
path lengths obtained for the WS and PN suggest some 
intrinsic property for each of those models as far as ran- 
dom paths are concerned. 

(e) Similarity of random path lengths in PA and ER 
models: This effect is particularly intense for m — 3 (see 
Figure [5^). Interestingly, these two models had been 
found [7[ to present similar values of several topologi- 
cal properties. Such an intrinsic structural congruence 



between the BA and PA models seems to have implied 
also in similar random path lengths. Given that the PA 
model is obtained by applying the star-path transforma- 
tion over a BA model, it is surprising that such a type 
of networks would result similar to the ER model in so 
many aspects. 

({)Large dispersion of random path lengths in the GG 
model: This effect is particularly interesting given that 
the way in which such networks are construct could be 
expected to result in nearly degree regular networks (i.e. 
most nodes with similar degrees because of the randomly 
uniform spatial distribution of nodes). However, the ran- 
dom paths obtained from the GG model exhibited widely 
varying lengths extending from too short (about 5 edges) 
to very long (over 80 edges). No other considered model 
yielded such a large variance of random path lengths. 
This property is interesting because it implies that it is 
in principle very difficult to foresee how well a geograph- 
ical model (at least of the GG type) can be effectively 



5 



covered by a moving agent. In order to understand bet- 
ter the reason for such a large variance of path lengths 
in this model, it would be interesting to consider con- 
figurations involving larger values of d (see Section . 



B. Convergence to the Maximum Random Path 



Having investigated random paths in the considered 
complex networks models, we now proceed to study the 
longest paths in those structures. 

Figure [3] illustrates the length of the longest random 
path obtained along the following of successive random 
paths in a single network of each considered type with 
TV = 100 and m = 3 considering the initial 200 ran- 
dom paths. It is clear from these examples that the 
length of the currently maximum random path tends to 
increase fast and steadily to a stable value, defining a 
plateau, which tends to correspond to the maximum ran- 
dom path length P. Experiments involving 2500 steps 
confirmed that the plateaux tend to remain stable, sug- 
gesting that the longest path has indeed been identified 
after 200 steps. Though the standard deviations vary 
from one case to another, they are mostly small (except 
for the GG networks). 

Figure m presents the longest path lengths for N = 100 
and m — 5 along the 200 steps (recall that each step 
involving obtaining a random path and comparing its 
length with that of the longest current path). Observe 
that the increase in m, implying higher average degree, 
tended to augment the length of the longest path in 
all cases, particularly for the GG structures where the 
longest path length almost doubled. It is also interesting 
to notice that the standard deviations tended to decrease 
in all models for this higher value of m. 

The above results suggest that the simple methodol- 
ogy adopted in this work for identification of the longest 
path tends to converge fast to the longest path (or at 
least to a path very similar to that in length) at least for 
the considered values of N and m. Also, the obtained 
results indicate that most models, except mainly for the 
BA networks, contain longest paths passing through the 
great majority of the network nodes. Several Hamilto- 
nian paths have been found for WS networks. 



C. Similarity between random paths and the 
maximum random path 

Though the method for longest path identification 
tends to converge fast and steadily to a putative longest 
path, possibly representing the longest one in the net- 
work, it is interesting to consider how much each of these 
putative longest paths followed during the algorithm re- 
sembles the longest path obtained for a very large number 



of steps (henceforth called the reference longest path). 
The importance of such an investigation becomes clear 
when we consider that a single network may have more 
than one longest path (with the same or very similar 
lengths), possibly even involving the same nodes. Pro- 
vided we can quantify the difference between the refer- 
ence longest path and each considered random path, the 
study of the variation of such a difference along the steps 
can supply a valuable information about the eventual 
multiplicity of longest paths in the network. As an exam- 
ple of a extreme situation, a completely connected net- 
work will contain a large number of Hamiltonian paths, 
i.e. paths covering all nodes in the network. 

In this section we report experiments aimed at quan- 
tifying the difference between the subsequent random 
paths and the reference longest path considering the over- 
lap of nodes and edges. 

Given two paths PI and P2, each with Nl and A^2 
nodes (and therefore El = A^l - 1 and E2 = N2 - 1 
edges) , it is possible to estimate the dissimilarity between 
these two paths through the following index 



" max{m,N2) ' ' 

where M is the number of nodes which are different 
in both paths. Therefore, the index D„ varies between 
and 1, achieving its minimum value for identical paths. 
However, such a minimum value does not imply that 
the two paths are identical (though they have the same 
nodes) , because each of them may involve different edges 
from the original network. 

It is also possible to define a dissimilarity index con- 
sidering the overlap of edges along the two paths, i.e. 



" max{El,E2) ^ ' 

where Q is the number of different edges in both paths. 
Again, we have that < 5*6 < 1. The minimum dissimi- 
larity between the edges of the two paths again provides 
no guarantee that the two paths are identical because 
the order of the edges is not taken into account in Equa- 
tion [2] Though it would be possible to obtain a com- 
pletely strict measurement of similarity between the two 
paths by considering distances between the two respec- 
tive adjacency matrices representing the paths, such a 
measurement would not work for cases where the two 
paths involve the same nodes but different edges. For 
such reasons, the similarities between each random path 
and the reference longest path are henceforth performed 
in terms of the dissimilarity indices given in Equations [1] 
and H 

Figure [5] shows the node dissimilarity for all cases, con- 
sidering iV = 100 and m = 3, along 500 steps (from step 
700 to 1200) of the algorithm for longest path estimation. 
The reference longest path corresponds to that obtained 
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FIG. 2: The average and standard deviations of the lengths of 2500 random paths selected for 50 realizations of each model 
considering A'' — 100 and m = 3. 



after 2500 steps in each case. The results in this figure in- 
dicate that each network model is characterized by differ- 
ent dismilarities. The smallest dissimilarities (about 7%) 
were obtained for the WS networks. Combined with the 
fact that the WS tends to lead to long random paths, the 
above result means that the longest paths in these struc- 
tures tend to incorporate the same nodes. The largest 
dissimilarities (both average and standard deviations) 
were obtained for the geographical model, reflecting the 
large dispersion of random path lengths already observed 
for that type of networks. Figure [5] shows the node dis- 
similarities obtained for to = 5, all of which smaller than 
the respective counterparts in Figure [5l 

The edge dissimilarities for n = 100 and m = 3 are 
shown in Figure [T] All such relative values are higher 
than those respectively obtained for the node dismilarity 
(please refer to figure O . It is also clear from Figure [7] 
that, except for the WS model, the edges dissimilarities 
were all close to 80%. Therefore, despite the fact that 
the random paths tend to share several nodes with the 
longest path in those networks, they involve significantly 
different edges. The edges similarities obtained for n — 
100 and to = 5 are presented in Figure [51 These results 
are remarkably similar to those in Figure [7l except for 
the lower dissimilarity obtained for the WS model. Such 
a result suggests that the average degree seems to affect 
relatively little the edge dissimilarities obtained for the 
considered models, at least for iV = 100. 

V. CONCLUDING REMARKS 

The node degree and shortest path length have re- 
ceived special attention from the complex network com- 
munity as measurements of the topology of complex net- 
works because of their intimate relationship with the crit- 



ical properties of scale-free and small-world, respectively. 
Though important, such measurements are not sufficient 
to describe the complete topology of complex networks, 
in the sense that many distinct networks may yield iden- 
tical average node degree or average shortest path length 
(e.g. Q). It is therefore important to devise different 
ways of looking at complex networks. While random 
walks have been largely applied in statistical physics and 
complex networks, random path- walks (random paths for 
short) have received substantially less attention (e.g. [8]). 
Paths are important because: (a) they necessarily ter- 
minate (a walk can proceed forever by passing through 
repeated edges and nodes); (b) paths are special cases 
of minimal spanning trees [l5|, therefore containing the 
smallest number of edges required to cover the respective 
nodes; and (c) they exist in smaller number than walks, 
as implied by the requirement of not repeating edges or 
nodes. The current work explored the possibility to in- 
vestigate the structural properties of complex networks 
as far as two types of paths are concerned: random and 
longest. 

Random paths are obtained by selecting a node ran- 
domly from the network and initiating a path from that 
node, until the path terminates at its both extremities. 
The statistical study of the distribution of type of paths 
is interesting because the variation of such structures can 
reflect intrinsic structural properties in the networks. In- 
deed, the experimental results reported and discussed in 
this article indicated that the several types of consid- 
ered networks — namely ER, BA, WS, GG, PA and PN 
— are characterized by distinctive average and standard 
deviation values of their random path lengths. The BA 
model, in particular, yielded the smallest random paths, 
which is explained by the fact that once a hub is visited 
it will block access to many other nodes to which it is 
connected. The longest random paths were obtained for 
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FIG. 3: The average ± standard deviation of the lengths of the longest path obtained for 50 realizations of each of the considered 
models with N = 100 and m = 3. Observe the steady convergence to the plateau associated to the longest path. 



the WS model. Random paths arising from PN networks 
tended to present remarkably similar lengths. The ran- 
dom path lengths tended to increased with the average 
degree in all cases. 

Largely overlooked in the literature, the longest path 
of a network constitutes a particularly interesting global 
underlying structure. Though the determination of the 
longest path in a network represents an NP-complete, we 
have shown in the current article at least for relatively 
small networks {N = 100) and considered average de- 
grees — that a simple algorithm involving the comparison 
of subsequent random paths tend to converge to a long 
path fast and steadily. Such an efficiency paves the way 
to comprehensive investigations of the distribution and 
properties of the longest paths in the several considered 
network models. Again, and for similar reasons, the BA 
structures tended to present the shortest longest paths, 
while the WS networks systematically yielded longest 
paths comprising all the nodes. The longest paths tended 



to increase with the average degree in all cases. At the 
same time, the longest path lengths were characterized 
by surprisingly small standard deviations, except for the 
GG model. 

Several are the future investigations motivated by the 
reported results. First, it would be interesting to in- 
vestigate the stochastic distribution of the random and 
longest paths considering alternative network models. It 
would be also interesting to try to relate intrinsic random 
and longest path properties to specific network measure- 
ments in addition to the average degree. In another di- 
rection, it would be relatively easy to extend the consid- 
ered concepts and methods to the analysis of random and 
longest cycles, which are closely related to random and 
longest paths. The concept of similarity between paths 
(as well as walks) can also be refined and explored further 
in order to quantify how diverse are the paths performed 
in different networks can be. Such investigations could 
provide valuable information about the transient dynam- 
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FIG. 4: The average ± standard deviation of the lengths of the longest path obtained for 50 realizations of each of the considered 
models with N = 100 and m = 5. Observe the effect of the increase in the average degree {{k) = 2m) in the extension of the 
longest path in all cases, particularly for the GG model. 



ics of the respective networks. 
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FIG. 7: The edge dissimilarity between the subsequent random paths at each step and the reference path for N = 100 and 
m = 5. 
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FIG. 8: The edge dissimilarity between the subsequent random paths at each step and the reference path for N = 100 and 
m = 5. 
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